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ABSTRACT 

We produce synthetic images and SEDs from radiation hydrodynamical simulations 
of radiatively driven implosion. The imaged bright rimmed clouds (BRCs) are mor- 
phologically similar to those actually observed in star forming regions. Using nebular 
diagnostic optical collisional line ratios, simulated Very Large Array (VLA) radio im- 
ages, Ha imaging and SED fitting we compute the neutral cloud and ionized boundary 
layer gas densities and temperatures and perform a virial stability analysis for each 
cloud. We determine that the neutral cloud temperatures derived by SED fitting are 
hotter than the dominant neutral cloud temperature by 1 — 2 K due to emission from 
warm dust. This translates into a change in the calculated cloud mass by 8 — 35 %. 
Using a constant mass conversion factor (C u ) for BRCs of different class is found to 
give rise to errors in the cloud mass of up to a factor of 3.6. The ionized boundary 
layer (IBL) electron temperature calculated using diagnostic line ratios is more accu- 
rate than assuming the canonical value adopted for radio diagnostics of 10 4 K. Both 
radio diagnostics and diagnostic line ratios are found to underestimate the electron 
density in the IBL. Each system is qualitatively correctly found to be in a state in 
which the pressure in the ionized boundary layer is greater than the supporting cloud 
pressure, implying that the objects are being compressed. We find that observationally 
derived mass loss estimates agree with those on the simulation grid and introduce the 
concept of using the mass loss flux to give an indication of the relative strength of 
photo-evaporative flow between clouds. The effect of beam size on these diagnostics 
in radio observations is found to be a mixing of the bright rim and ambient cloud and 
HII region fluxes, which leads to an underestimate of the cloud properties relative to 
a control diagnostic. 

Key words: stars: formation - ISM: HII regions - ISM: kinematics and dynamics - 
ISM: clouds - methods: observational - methods: numerical 



1 INTRODUCTION 

A central problem in the study of star formation and galaxy 
evolution is the effect of stellar feedback on star formation 
rates, efficiencies, and the initial mass function (IMF) (e.g. 
Dobbs et al.||2011[ |Elmegreenl[20TT| [Batel[20T2l |Elmegreen| 



2012 



Kennicutt fc Evans|2012[ and references therein) . Stel- 



lar winds drive into the surrounding material which is ex- 
pected to either induce star formation by triggering collapse 
in otherwise stable conglomerations or to hinder star for- 
mation by dispersing material and driving turbulence (ex- 
ample theoretical, observational and synthetic observational 
2001) JSmith fc Rosen||2005l jArce] 



studies include Matzner 



et al. 2010 Offner et al 



2011). Supernova feedback is an- 



other kinematic mechanism whereby the explosions of mas- 
sive stars locally disperse potential star forming material, 
but may cause large-scale collapse of nearby gas clouds (see 
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Hensler 2011 and references therein). Young stars, specif- 
ically the most massive (i.e. OB), emit large amounts of 
ionizing radiation into the surrounding gas. This quickly 
ionizes and heats the surroundings, generating shocks that 
sweep up material into a shell that can either become locally 
gravitationally unstable (the collect and collapse scenario - 
|Deharveng et al.|2005||Zavagno et al.|2006||Dale et al.|2007| ) 
or drive into pre-existing density structures and compress 
them to form stars, this final process is known as radiatively 
driven implosion (RDI). 

RDI has been subject to a large amount of numerical 
modelling which has repeatedly demonstrated that the com- 
pression and implosion of clouds is possible (e.g 
|et al.p982l|Bertoldi|1989||Lefloch 
Peynet fc Burkert|2003||Miao et al 



Sandford 



Lazarcff 1994| |Kessel- 
2009[|Gritschneder et al 
Bisbas et al.||2011| |Dale fc Bonnell||2012[ |Haworth fc 



Harries|2012 1. The resulting objects are usually bow shaped 



thin dense shells or cometary structures, in qualitative agree- 
ment with studies such as the northern hemisphere bright 
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rimmed cloud (BRC) survey of Sugitani et al. ( 1991 1, which 



found that BRCs can be broadly classified by their rim mor- 
phology, being either: type A (moderately curved), type B 
(tightly curved) or type C (cometary). 

A number of bright rimmed clouds have been identi- 
fied in the vicinity of which are often collections of young 
stars (studies include Sugitani et al.|1991||Ogura et al.|2002 



Thompson et al. 200 4||Lee fc Chen|2007||Beltran et al. 2009 
Morg an et al.||20 09 Chauhan ct al. 2009[ |Choudlmry et al.| 



2010 Chauhan et al. 2011 1, these are often cited as examples 



of triggered star formation via RDI and have been studied 
using a range of techniques. One of the most popular signa- 
tures of triggered star formation is an age gradient in YSOs, 
with bluer and therefore supposedly older sources located in 
closer proximity to the triggering star. |Beltran et al.| ( |2009[ ) 
study the BRC IC 1396N in the Cep OB2 association and 
caution against this approach, as a collection of stars of sim- 
ilar age may exhibit an apparent observational age gradient 
if they are sequentially exposed to high ionizing flux and 
stripped of their circumstellar environment. 

An alternative to studying stellar age gradients is to an- 
alyze the neutral cloud and ambient ionized gas properties. 
|Lefloch fc Laza rcff (1994) calculated semi-analytic two di- 
mensional models of RDI and included figures of the system 
emissivity which were found to have similar morphology to 
the BRC types observed in star forming regions. They fur- 
ther compared the cloud and ionized boundary layer (IBL) 
pressures using a multiwavelength molecular line and ra- 
dio emission study of CG5 in IC1848 ( |Lefloch et aT1|l997| | 
and found greater pressures in the IBL by approximately a 
factor of 10, indicative of shock compression. They also cal- 
culated an estimate of the incident ionizing radiation flux 
and mass loss rates, finding values of 4.8 x 10 9 cm -2 s -1 



and 105 M Myr -1 respectively. Thompson et al. (20041 and 
|Morgan et al.| ( |2004[ ) performed a similar analysis on a larger 
number of clouds, finding that for only a small number the 
IBL is at higher pressure and that more frequently the ob- 
jects are in pressure balance. The aforementioned studies 
suggest that greater external pressure will lead to the even- 
tual dispersal of the cloud, whereas clouds that are overpres- 
sured relative to the surroundings will stall the shock until 
enough material is accumulated and the pressure sufficiently 
increased to continue driving into the cloud. The presence 
of photo-evaporative flow in Lefloch & Lazareff ( 1994 1 and 



Thompson et al. ( 2004| is also an important factor in iden- 
tifying RDI. In addition to calculating the mass loss rate, 
striations in Ha imaging about the BRC is identified as an 
indication of such a flow. As yet, there is no direct transla- 
tion between numerical models and these observational di- 
agnostics of BRCs. 

In this paper we use the final state of the RDI models 
generated in |Haworth fc Harries| ( |2012| ) that incorporated 
the diffuse field. These calculations considered the RDI of a 
Bonnor-Ebert sphere (BES) at three distances from the trig- 
gering star and gave rise to clouds that appear to be type 
A (the low flux model), B (the high flux model) and type 
B-C (the medium flux model). These models also vary in 
the nature of the driving shock, with the medium and low 
flux models being driven by a photoionizing shock/photo- 
evaporative flow and the high flux model in pressure balance. 
With these final grid states we perform synthetic Ha and ra- 
dio continuum imaging, as well as generate synthetic SEDs 



to perform standard diagnostics that look for RDI. We also 
explore the use of long slit spectroscopy to calculate diagnos- 
tic forbidden line ratios, determine the nebular conditions 
and infer the stability of the BRCs. In addition we investi- 
gate the observational characteristics of photo-evaporative 
flows and mass loss rate calculations. This is all performed 
from the unique perspective in which the actual conditions 
of the system are known. Through this process we aim to 
guide observers towards unambiguous signatures of triggered 
star formation in BRCs and to investigate the accuracy and 
applicability of these standard diagnostic techniques. 



2 NUMERICAL METHOD 
2.1 Overview 

We used the grid based radiation transport and hydrody- 



namics code torus (Harries 2000 Acreman et al. 2010 



Harries|[2011 Haworth & Harries 2012) to perform the cal- 
culations in this paper. Using the final states of the mod- 
els which included the diffuse field in Haworth & Harries 



(20121 we perform photoionization calculations with addi- 



tional atomic species. Dust is then added to the resulting 
grids as part of post processing and the final grid states are 
used to generate SEDs and images for diagnostics. 



The models from Haworth & Harries (20121 comprised 



a star at varying distances from a BES. Radiation from the 
star (which was off the edge of the grid) entered the grid at 
the —x boundary in a plane parallel manner, under the as- 
sumption that the grid was sufficiently far from the ionizing 
star for this to be the case. An ionization front formed and 
accumulated material that was driven into the BES, com- 
pressing it. Material was allowed to stream freely off the 
edge of the ±x boundaries but not allowed to re-enter the 
grid. The other boundaries were periodic. 



2.2 Photoionization 

The photoionization routine follows that detailed in |Er- 



colano et al. ( 2003 1, Wood et al. ( 2004 I and Haworth & Har- 



(20121 and involves propagating photon packets from 



sources through a computational grid. A photon packet is a 
collection of photons for which the total energy is constant, 
being the total luminosity of the stellar photon sources di- 
vided by the total number of packets. The number of pho- 
tons in the packet varies with frequency. As photon packets 
traverse cells on the grid they modify the energy density 
in the cell region, allowing the ionization and temperature 
states to be calculated. The diffuse field is included in mod- 
els in this paper, whereby a new photon packet is immedi- 
ately generated upon absorption with a new frequency and 
direction and is further propagated over the computational 
domain. This process repeats until the packet escapes the 
grid. We also include periodic photon packet boundary con- 
ditions, without which the ionizing flux at the edges of the 
grid can be underestimated if the model geometry is peri- 
odic. However, in order to avoid long loops for low frequency 
packets which have a small probability of interaction, each 
photon packet is only allowed to cross a periodic boundary 
once. Ionization balance in all calculations is determined by 



solving the ionization balance equation ( Osterbrock 1989 1 
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and photoionization equilibrium is assumed. In this paper 
we include treatment of a range of atomic constituents: hy- 
drogen, helium, carbon, nitrogen, oxygen, neon and sulphur. 
Thermal balance is calculated by finding the temperature at 
which the heating and cooling rates match in the manner de- 



scribed in Wood et al. (2004 1. Hydrogen and helium heating 
is treated as well as recombination, free-free and collisional 
line cooling. 



2.3 Synthetic Imaging and SEDs 

Synthetic images are generated by accumulating photon 
packets from the grid in a 2D array of collecting bins. The 
basic method follows the Monte Carlo scheme detailed in 



Hillier (19911 and Harries (2000|, an example application is 
Kurosawa et al.| ( |2004[ ). For an observer location (which can 
be specified arbitrarily) relative to the grid, a pixel array 
is generated. The emissivity across the grid is then calcu- 
lated and can include dust continuum, free-free continuum, 
forbidden and recombination line radiation. The image is 
constructed by propagation of photon packets from dust, 
gaseous and stellar sources, with each packet carrying power 
P given by 



P = 



U + W t ot 

■^V montc 



(i) 



where L* , Wtot and iV mc ,ntc are the total stellar source emis- 
sion, the total emission from the gas and dust and the num- 
ber of photon packets used in the image generation respec- 
tively. 

The Monte Carlo method is subject to Poisson noise 
analogous to that inherent in the collection of photons in real 
imaging, therefore we employ variance reduction techniques 
to reduce the calculation time required to reach adequate 
signal to noise. To improve sampling in the gas the proba- 
bility of a photon packet originating from a stellar source is 
fixed at 0.1 and the packet assigned an appropriate weight 
w. Following initial generation and each scattering event of 
photon packets, an additional packet is forcibly directed to- 
wards the observer with a modified weighting the, so called, 



'peel-off' technique (Cashwell & Everett||1959 Yusef-Zadeh 



et al. 19841. Once a photon packet escapes the grid in the 



direction of the observer it changes the flux F 7 in whichever 
pixel it intersects to 



F- 



7 + 1 



Pwe 
An 



(2) 



where w is the photon packet weight and r is the optical 
depth along the packet's path to the observer. The accumu- 
lated photon packet contributions are converted into a final 
image in distance independent units of mega Janskys per 
steradian. 

Spectral energy distributions are also generated via 
Monte Carlo radiative transfer in a similar manner. Rather 
than contributing to a 2D pixel array photon packet intensi- 



ties are binned by frequency to form the SED (see Kurosawa 
et al.|2004| >. 



2. 3. 1 Dust properties 

Dust is not directly included in the photoionization calcula- 
tions in this paper, rather it is added to all cells where the 
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Figure 1. The dust opacity as a function of wavelength. 



temperature is less than 1500 K prior to image/SED gener- 



ation (e.g. Dullemond et al. 20011. Unless otherwise spec- 
ified, a dust to gas ratio of 1 x 10 -2 is used in all images 



generated in this paper. We assume spherical silicate dust 
particles that follow a standard interstellar medium power- 
law size distribution (e.g. Mathis et al. 19771. The optical 



Draine fc Lee| (|1984). We use a 



constants are taken from 

pre-tabulated Mie-scattering phase matrix. The wavelength 
dependency of the dust opacity is given in Figure[l] Dust was 
not included in the original radiation hydrodynamic models 
because a simplified thermal balance calculation was used 
to remain consistent with previous models. A realistic treat- 
ment of the dust would therefore not be achieved. It is ex- 
pected that dust would not significantly modify the radia- 
tion hydrodynamic calculations as it would be destroyed in 
the ionized medium, not altering the pressure external to 
the cloud and not increasing the supporting cloud pressure. 
However confirming this is a subject for future study. 



2.4 Radio Imaging 

Synthetic images generated at radio wavelengths using the 
method described in section [273] will have a sub-pixel beam 
and a Poisson noise level determined by the number of pho- 
ton packets used in the image generation, at a level that is 
possibly lower than the Gaussian noise associated with radio 
imaging. We therefore modify images at these wavelengths 
as part of post processing. The image is smoothed to a Gaus- 



sian beam using ACONVOLVE from CIAO v4.1 (Fruscione et al. 
|2006[ ) to a size appropriate to the half power beamwidth 
(HPBW) of the simulated instrument. We assume that the 
Gaussian smoothing outlined above extinguishes the effects 
of Poisson noise (which is relatively small) and then add 
Gaussian noise to the image separately using the STARLINK 
ADDNOISE routine. The emulated exposure time, beam size 
and noise level are chosen to be typical of those for a selected 
instrument. 
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3 DIAGNOSTIC TECHNIQUES 

To determine what an observer would infer following obser- 
vations of our model cloud we employ a range of standard 
diagnostic techniques that are used to determine the ionized 
and neutral gas properties. 



assume to be 1.36 in neutral gas following, e.g. |Hildebrand| 
(1983[);|Blitz fc Rosolowsky] (|2006 1. 



The model clouds are not supported by turbulent mo- 
tions owing to the ideal starting conditions. We assume that 
the sound speed in the neutral gas is isothermal and is cal- 
culated using the fitted dust temperature. 



3.1 Neutral cloud properties from SED fits 

In this work no molecular lines are treated. We therefore 
obtain the neutral cloud gas properties by fitting the sub- 
millimetre thermal dust continuum tail of the object SED 
with a greybody function of the form 



F„=nB„(T d ) (l-e" r -) 



(3) 



where Q, B v (T d ) and t„ are the solid angle subtended by the 
region for which the temperature is being derived, the fre- 
quency specific Planck function at dust temperature Td and 



the optical depth at frequency v respectively (Hildebrand 



1983 



|Dent et al.|1998||Thompson et al.|2004| [Morgan et al 



20081. We adopt a procedure in which the optical depth at 
a given frequency is paramaterised in terms of a reference 
frequency and optical depth i.e. 



— 7~ref 



^rcf 



P 



(4) 



(e.g. |Hildebrand| 1 1983] ) . We choose a reference wavelength 
of 850 /im and adopt a value of 2 for j3, the index specify- 
ing the frequency dependency of the dust emissivity, follow- 
ing |Hildebrand| ( p83) , [Thorn etal] pOOlf and |Morgan] 



et al. (20081. The reference optical depth is calculated us- 



ing the equation for submillimetre optical depth from |Hilde-| 
brand| ( [T^83| > 



F v WelB v {T d ) 



(5) 



where #r is the angular radius of the region over which the 
flux is being integrated in radians. We fit the SED using 
a chi-square minimization. Although we can fit across an 
entire calculated spectrum, we typically do so over only a 
small range in the manner of observational studies such as 



Thompson et al. ( 2004 1. The cloud mass is then found using 



the established method of Hildebrand (19831, whereby the 



total gas and dust mass of the cloud is given by 



M = 



d 2 F v C„ 
B„(T d ) 



(6) 



where d is the distance of the cloud from the observer and 
C v is a mass conversion factor. A value for C v needs to 



be appropriately selected (Hildebrand 1983 Draine & Lee 



1984| |Ossenkopf fc Henning|1994||Kerton et al.|2001[ ). Typ 

ically, smaller values of C v apply to cold, high density, re- 
gions and larger values to low density regions (e.g. the diffuse 



ISM). As an example, the value selected by Thompson et al. 
( |2004[ ) for 850 /im analysis of clouds with assumed densities 
of 10" cm" 3 is 50 gem -2 . An expression for C„ is given by 



Hildebrand| ( |1983[ ) 



CV = [Nh/t v ] mH/i 



(7) 



where Nn is the column density, ttih the hydrogen mass and 
H is the mean particle mass relative to hydrogen, which we 



3.2 Ionizing flux, mass loss rates and IBL electron 
densities from radio emission 

We calculate the conditions in the IBL of the BRC following 



the standard radio diagnostics of Lefloch et al. ( 1997 1 which 



are used in, for example, Morgan et al. (20041, Thompson 



et al. (20041 and Urquhart et al. (2006). This technique does 



not constrain the temperature in the IBL so we assume the 
standard value of 10 4 K following the aforementioned stud- 
ies. 

An observational estimate of the ionizing flux $ per 
square centimetre per second impinging upon the BRC can 
be made from the 20 cm free-free emission integrated flux 
from the IBL by rearranging equation 6 from |Lefloch et aL~| 
(T997] 



$ = 1.24 x 10 10 P„T o ' 35 j/ 01 * 



(8) 



where T e is the electron temperature of the ionized gas in 
K, 6 is the angular diameter in arcseconds of the region over 
which the flux is integrated and F v is the integrated flux in 
mjy at frequency v in GHz. 

The electron density can also be calculated using 



n c = 122.41 



^7x1.35^0.151- 



t]R 



(9) 



where r\ is the thickness of the ionized boundary layer as 
a fraction of the cloud radius and R is the cloud radius in 
parsecs ( |Lefloch et al.|1997[ ). [Bertoldi|||1989|| give r\ typically 
« 0.2 which is the value adopted by Morgan et al. (20041, 
|Thompson et al.| ( |2T)0"4"| ) and |Urquhart et al.| ( |2006[ ), we there- 
fore also assume this value. The electron density allows the 
ionized boundary layer gas pressure to be derived using 



Pi = 2p ; c 2 



(10) 



where Pi, pi and c\ are the pressure, density and sound 



speed in the ionized boundary layer (e.g. Morgan et al 



2004 Thompson et al. 20041. We assume that in the ion- 



ized regions where these diagnostics are applied that pi = 
n e mn and that the sound speed is isothermal, i.e. ci = 
\Jk-B.T c / [im-n, where ks is the Boltzmann constant and \x is 
the mean particle mass relative to hydrogen, here assumed 
to be 0.6 in the IBL, the value for ionized gas of solar com- 
position. 

The mass loss rate from the cloud due to ionized gas 
streaming away from the surface (photo-evaporative flow) is 
calculated using 



M = 4.4 x 10~ 3 $ 1/2 P 3/2 M Q Myr 



(11) 



where R is the cloud radius in parsecs (Lefloch & Lazareff 
|1994| [Thomps on et al.|2004| ) . The strong photo-evaporative 
flows identified in some of the RDI model clouds from Ha- 



worth & Harries (2012), which are also studied here, will 



provide a valuable test of the mass-loss rate estimate from 
Equation |11| 
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3.3 H a emission 



We generate images of H a recombination emission as it has 
previously been used to identify photo-evaporative flow via 
striations perpendicular to the bright rim ( jThompson et al~ 
20041. Imaging at this wavelength (6563 A) is also useful for 



tracing ionized gas and is regularly used as a basis for iden- 
tifying regions over which to perform long slit spectroscopy 
(e.g. |Tiillmann et al.||2003[ ) . 



3.4 



IBL properties from optical collisional line 
ratios 



The conditions of HII regions and planetary nebulae have 
long been studied using ratios of forbidden lines that are 
sensitive to the electron density or temperature (e.g. |Os- 



terbrock|1989| |Spitzer||1998| |Caplan et al.pOOOl |Deharveng 
et al.|2000 Lagrois et al.|2012| and references therein). Sin- 
gle forbidden lines have been used in studies to identify pre 
main sequence (PMS) stars in and around BRCs as an in- 
dicator of triggered star formation (e.g. |Lee et al.|2005| |Lee 
& Chen 20071. However, the line ratios have not yet been 



applied to calculating the IBL properties and hence the rela- 
tive pressures of the IBL and the neutral cloud. We therefore 
derive the ionized gas conditions using these diagnostic ra- 
tios, in addition to the radio analysis, to assess their future 
use as a tool for identifying RDI. 

Forbidden line intensities are usually obtained by per- 



forming slit spectroscopy of the system (e.g. Tiillmann et al. 
|2003[ ). We emulate this form of analysis by generating syn- 
thetic images at a range of forbidden line wavelengths and 
calculating the intensity across a pseudo-slit region on the 
image. This is akin to choosing the intensity at a specific 
wavelength on the slit spectrum. 

We calculate the electron density using the [O II] 
3729A/3726A ratio by linearly interpolating between the 
tabulated ratio and electron density values given in Table 2 
of Wang et al.| (|2004[ > as well as some values from Table 5.2 
of Osterbrock (119891). We assume a maximum ratio of 1.5, 



determined by the ratio of the level statistical weights ( Os 



terbrock||l989 1 and use a corresponding minimum electron 
density of 10 cm -3 . The resulting variation in the collisonal 
line ratio as a function of logarithmic electron density is 
given in Figure [2] This electron density value is then used 
in the [O III] line ratio 



(12) 



JA4959 + JA5007 _ 


7.90exp(3.29 


x 10 4 /T) 


,?A4363 


1 + 4.5 x 10" 




the [N II] line ratio 






JA6548 + JA6583 


8.23exp(2.50 


x 10 4 /T) 


JA5755 


1 + 4.4 x 10" 


3 n e /TV2 


and the [Ne III] line ratio 






JA3869 + JA3968 _ 


13.7exp(4.30 


x 10 4 /T) 



JA3343 



(13) 



(14) 

1 + 3.8 x lO-^c/TVa v ; 

expressions to find the 



low density limit (riH < 10 cm ) 
temperature (e.g. |Osterbrock| 1989"] ) 

Hence from one electron density sensitive line ratio and 
one temperature sensitive ratio an estimate of the electron 
density and temperature in the ionized gas can be obtained. 
With an estimate of the temperature and electron densities 



co 

<M 
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Figure 2. The [O II] (3729 A/3726 A) ratio variation with loga- 
rithmic electron density. 



the pressure is calculated using equation |10| This method 
has the advantage that the temperature of the ionized gas 
is determined directly, rather than being assumed to be 



10000 K as is the case for the radio diagnostics (e.g. Thomp- 
son et al.||2004[ ). 



4 TESTING 

The TORUS radiation hydrodynamics scheme is tested exten- 



sively in Haworth & Harries (20121. There are also numer- 



ous applications of the imaging and SED generating routines 
available in the literature (e.g. Harries et al.|200~ Kurosawa 
et al.(2004 i. Tests of the ratio diagnostics and the SED fit- 
ting are given here. 



4.1 Diagnostic line ratio testing: the HII40 
Lexington benchmark 

The HII40 Lexington benchmark involves, in one dimension 
and assuming no time evolution, calculating the tempera- 
ture and ionization structure of the constant density gas 
surrounding a star at 40000 K ( |Ferland|[l995l ). The test in- 
cludes no dust and the geometry of the cloud is a thick, 
constant density, spherical shell centered on the star with 
an inner radius 0.97 pc (3 x 10 18 cm) from the star. A list 
of the parameters used in the HII40 benchmark is given in 
Table [l] Given the well studied nature of this temperature 
distri bution ( |Ferland|1995"||Ercolano et al.|2003"||Wood et al. 
2004[ ) it is an ideal test of the accuracy of our forbidden line 
ratio diagnostics. We perform a photoionization calculation 
of a three dimensional version of the HII40 benchmark with 
the ionizing star set at the grid centre and apply forbidden 
line ratio diagnostics to synthetic images of the converged 
system. Given the large amount of neutral foreground ma- 
terial and the low dust fraction in the hot HII region, the 
dust to gas ratio in this test is set to a negligibly small value 
so that extinction does not have to be accounted for. The 
model resolution, number of image pixels and slit size are 
chosen to match that of the RDI grids, images and emulated 
slit spectroscopy of BRCs used later in this paper. The pa- 
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Table 1. Parameters for the HII40 Lexington benchmark which 
is used to test diagnostic line ratios. 



Variable (Unit) 


Value 


Description 


A eft 


40000 


< ~Jr»iiTT , o ^*fT^ l ^'i~^^r^* t"OTYii~ir>T"ntiiT'f"> 

kJUIilLL LllljL LI VL LLjlll IJ y^L CX\i Lll L 


it* ^1X0 ) 


18.67 


Source radius 


n H (cm" 3 ) 


100 


Hydrogen number density 


logi (He/H) 


-1 


Helium abundance 


logi (C/H) 


-3.66 


Carbon abundance 


logi (N/H) 


-4.40 


Nitrogen abundance 


logi (O/H) 


-3.48 


Oxygen abundance 


logi (Ne/H) 


-4.30 


Neon abundance 


logi (S/H) 


-5.05 


Sulphur abundance 


L (pc 3 ) 


16 3 


Computational domain volume 


^cclls 


128 3 


Number of grid cells 




401 2 


Synthetic image pixels 




10 8 


Photon packets used in 






synthetic image generation 


L B (pixels) 


100 


Spectroscopic slit length 


W B (pixels) 


2 


Spectroscopic slit width 



rameters associated with synthetic imaging of this test are 
also given in Figure [l] 

The converged temperature state of the grid is given in 
Figure [3] and a colour composite image comprising Ha (red) 
the 5007 A [O III] line (green) and the 3968 A [Ne III] line 
(blue) is given in Figure [4] upon which is marked the region 
covered by the spectroscopic slit. In addition to the slit loca- 
tion in Figure [4] we repeat the diagnostics for slit positions 
at ±1,2 pixels in the a;-direction to determine the level of 
uncertainty when using a single slit position to infer the con- 
ditions. The average electron density calculated using the [O 

II] diagnostic ratio and temperature calculated using the [O 

III] , [N II] and [Ne III] diagnostic ratios are given in Ta- 
ble[2] The prescribed hydrogen number density is 100 cm -3 , 
which is slightly higher than the inferred electron density. 
The HI fraction throughout the grid is not uniformly zero, 
however this accounts for less than 1cm -3 . Comparing the 
linear interpolation between our [O II] data points with a 
more sophisticated fit also only leads to an improvement of 
3 cm -3 . The remaining discrepancy is therefore believed to 
be characteristic of the diagnostic ratio. 

Although there is some variation between the calculated 
temperatures, they all lie within the 7000-10000 K range ex- 
pected from Figure [3] The actual average HII region tem- 
perature in the region covered by the slit is 7769 K, within 
600 K of each diagnostic and within 85 K of the diagnostic 
average. Based on the uncertainties in individual diagnostics 
calculated in this test, further electron densities derived in 
this paper using the [O II] line ratio are conservatively cal- 
culated to the nearest 10 cm -3 and individual temperatures 
derived using the [O III], [N II] and [Ne III] line ratios are 
calculated to the nearest 10 K. 

4.2 SED greybody fit testing: a cold uniform 
sphere 

To test the cloud masses and temperatures derived by grey- 
body SED fitting we consider the simple system that is a 
three dimensional uniform density sphere at constant tem- 
perature, surrounded by vacuum. We choose a density of 
100 hydrogen atoms per cubic centimetre and a cloud ra- 



Table 2. The Lexington HII40 test average conditions from di- 
agnostics using a spectroscopic slit at 5 locations. 

Ratio n c (cm -3 ) T e (K) 

[OH],i3729/j3726 88 ± 2 

[O III], 05007 + j4959)/i4363 - 7215±6 

[Nil], 6583 +i6548)/i5755 - 8366±| 

[Ne III] , (j 3 g 68 + J 3 869)/i3343 - 7474±| 




Figure 4. A colour composite image of the HII40 test using Ha 
(red), the [O III] 5007 A line (green) and the [Ne III] 3968 A line 
(blue). Overlaid is the slit region used to emulate slit spectroscopy. 
The image side spans 16 pc 
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Figure 5. The uniform cloud model SED (red line) and the grey- 
body fit (blue crosses). 



Table 3. Cold uniform sphere fitted masses. 



d, (gem" 2 ) 


Mass M 


50 


9.9 


100 


19.8 


200 


39.6 


214 


42.4 


250 


49.5 



dius of 1.6 pc (equivalent to the cut-off radius of the start- 
ing Bonnor-Ebert sphere in the models of Hawort h fc Harries] 
2012[ ). The surrounding vacuum is prescribed a number den- 
sity of 10 -30 hydrogen atoms per cubic centimetre so that 
it will only negligibly contribute to the system SED. The 
model is prescribed a uniform temperature distribution of 
10 K throughout the grid. 

We produce an SED for an observer situated 1000 pc 
from the grid. The spectrum is fitted in the manner de- 
scribed in section |3.1| from 450 fim to 850 /im following 
Thompso n et al.] ( |2004[ ) and is shown in Figure [5] The re- 
sulting calculated temperature is 9.5 K, within 5% of the 
prescribed value. The total mass of this cloud is 42.4 Mq. 
We use a range of values for C„ to calculate a fitted mass 
from the SED, the results of which are given in Table [3] We 
find that the most appropriate value of C„ for number den- 
sities of order 100cm -3 at low temperatures is 214gcm~ 2 . 
The average column density, calculated by integrating the 
column density over the uniform sphere and dividing by the 
projected area, is 1.3 x 10 20 cm -2 , the optical depth calcu- 
lated using equation [H] is 2.2 x 10~ 6 and the corresponding 
C u value using equation [7] is 131 gem -2 , just under a fac- 
tor of 2 smaller than our inferred value. The reason for the 
discrepancy between the theoretical and inferred values is 
likely due to averaging the large variations in column den- 
sity across the sphere, from its maximum value at the cen- 
tre to approximately zero at the edge. Given that densities 
of order lOOmHcm" 3 are typical of the neutral clouds in 
the RDI models, we retain the presently inferred value of 
C v — 214 gem -2 in further SED fitting in this paper. 



5 RESULTS AND DISCUSSION 
5.1 The numerical models 

We add atomic species to the final state (after 200 kyr of 
evolution) of the RDI models in |Haworth fc Harries| ( |2012[ ) 
and perform photoionization calculations. The three models, 
each of which comprised a Bonnor-Ebert sphere (BES) ex- 
posed to a different level of ionizing flux (high, medium and 
low), underwent thermal and photo-evaporative compres- 
sion. The low flux model exhibited thermal compression of 
the BES and intermediate strength photo-evaporative flow, 
resulting in a type A bow. Instabilities in the ionization front 
also led to the formation of finger like objects in the wings 
of the low flux model. The medium flux model quickly ac- 
cumulated a dense shell of material resulting in a strong 
photo-evaporative flow that compressed the cloud to a type 
B-C cometary bow. The high flux model rapidly established 
pressure balance and exhibited only weak photo-evaporative 
flow, leaving what resembles a type B bow. A summary of 
the key parameters of these models is given in Table [4] 

The temperature state of each model following the pho- 
toionization calculation with additional atomic chemistry is 
given in Figure [6] In each model, typical ionized gas temper- 
atures range from 7000 — 10000 K and neutral temperatures 
are essentially uniform at 10 K. Of the cells with an HI ion- 
ization fraction greater than 0.5 and temperature less than 
1500 K: 97%, 96% and 98% are at 10 K for the low, medium 
and high flux models respectively. By mass 96%, 95% and 
98% of the cloud is at 10 K for the low, medium and high flux 
models respectively. Hot spots in the neutral gas of the bot- 
tom frame of Figure [6] are artifacts due to Monte Carlo sam- 
pling where, on rare occasion, high energy photon packets 
propagate a significant distance into the neutral gas. How- 
ever, these spots are too hot for dust to survive and are away 
from the region over which we analyze the IBL so will not 
affect our synthetic SEDs or long-slit spectroscopy. 

The low and medium flux models have retained a cold 
neutral gas morphology similar to that at the end of their ra- 
diation hydrodynamic evolution. Because the disruption to 
the original BES significantly modified the density distribu- 
tion in the low and medium flux regimes by excavating and 
accumulating material, the addition of metals (and there- 
fore forbidden line cooling) and more sophisticated thermal 
balance have little effect on the extent of the ionized region. 
This is not the case for the high flux model, which weakly 
modified the density distribution, achieving early pressure 
balance and establishing only weak photo-evaporative flow. 
As a result the enhanced cooling has shifted the ionization 
front away from the cloud structure that was formed in the 
radiation hydrodynamics calculation to a new position sim- 
ilar to that of the low flux model. This additional cooling in 
weakly disrupted gas is also responsible for the cool regions 
towards the edge of the medium flux grid. These alterations 
to the extent of the ionized regions hint at the importance of 
including additional atomic processes such as forbidden line 
cooling in radiation hydrodynamics calculations. Although 
the effect is small in the low and medium flux model, any 
diagnostic procedure applied to the high flux model will not 
be examining the pressure balanced system generated in |Ha-| 
|worth &: Harries] ( |2012[ ), but a new one in which D-type 
expansion about the HII region is yet to occur. 

The low and medium flux models both exhibit slightly 
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cooler regions where strong photo-evaporative flow occurs. 
This is the yellow region mirroring the central part of the 
bow of the top frame of Figure [6] for the low flux model 
and the shovel shaped yellow region about the tip of the 
cometary object in the middle frame of Figure [6] for the 
medium flux model. 



5.2 Optical Image Morphology 

A synthetic image of the low flux model using Ha (red), the 
5007 A [O III] line (green) and the 3968 A [Ne III] line (blue) 
is given in the top frame of Figure [7| The object clearly re- 
sembles a class A, moderately curved, bright rimmed cloud. 
Ha dominates the emission along the rim, making it ap- 
pear red. The signature of the photo-evaporative flow is the 
slightly darker region opposite the bright rim, which has 
been excavated to slightly lower densities due to the motion 
of the flow. The hot (blue) spots are due to Monte Carlo 
sampling as explained in section [5~T] From this edge-on in- 
clination the fingers due to instabilities (see section 5.1 1 are 



not clear because along a given line of sight there are multi- 
ple trunks and the average column density is approximately 
constant. Furthermore, the enhanced cooling due to forbid- 
den line processes means that the ionized regions between 
the trunks are now neutral, suppressing emission from them. 
A two colour image in the electron density sensitive lines of 
[O II] where the observer is inclined by 30 degrees relative 
to the ionization front is shown in Figure [8] Here the fin- 
gers and dark photo-evaporative flow region are much more 
clearly illustrated. The morphological resemblance between 
this image and that in Figure 1 of |Chauhan et al.| $20lT\ 
is notable, despite the fact that the additional cooling has 



made the fingers harder to discern. In Chauhan et al. (2011 1 



it was suggested that the finger like objects in the wings of 
the BRC were due to instability, this image and the radia- 
tion hydrodynamic model it is generated from qualitatively 
support this hypothesis. 

A synthetic image of the medium flux model using Ha 
(red), the 5007 A [O III] line (green) and the 3968 A [Ne III] 
line (blue) is given in the middle frame of Figure [7] This ob- 
ject resembles a type B-C BRC, having an almost cometary 
appearance that is again dominated by Ha emission. The 
ribbing effect interior to the cometary structure is a result of 
the low resolution used in the radiation hydrodynamics cal- 
culation and has no physical significance. The bright ribs are 
the cells exposed to the direct stellar radiation field and the 
darker regions are shielded cells, for which radiative heating 
comes primarily from the diffuse field. Despite being known 
to have a stronger photo-evaporative flow than the low flux 
model, it is difficult to distinguish between a dark excavated 
region and neutral foreground material. Again, we take an 
image at thirty degrees inclined relative to the ionization 
front to illustrate this point in Figure |9j in which the dark 
foreground material starts moving out of view. This image 
suggests that the BRC is situated in the bottom of a basin of 
ionized gas. Furthermore, a darker region excavated by the 
photo-evaporative flow is visible around the BRC in these 
electron density sensitive diagnostic lines. 

A synthetic image of the high flux model using Ha (red), 
the 5007 A [O III] line (green) and the 3968 A [Ne III] line 
(blue) is given in the bottom frame of Figure [7] Note that 
the radiation hydrodynamics model did not significantly al- 




Figure 6. The temperature distribution (in Kelvin) for each 
model, low flux to high flux from top to bottom. The grid is 
4.87 pc to a side. 



Triggered star formation diagnostics 9 



Table 4. Key RDI model parameters. 



Variable (Unit) 


Value 


Description 


Di (r>c) 
-^low vPW 


(—in fis n 0) 


Source position (low flux) 


^mcd \PW 


(—A 7s n n\ 


Sniivro nfwi ("i/in 1 tti r>rl 1 1 1 ill fluvl 

. _J l J 1 1 1 lA UUS1L1U11 l 111(^111 Hill 11 UX 1 


-^hi \t ,Kj J 


f_a "?8 0") 


Source position (high flux) 




40000 


. H 1 1 1 1 ' ' L11LLL1V1 LL111 1H_1 d L 111 




10 


"^niiT'fT 1 ruiliiiQ 


<&i ( cm ~ ^ c; - 1 ^ 
^low V*" 111 s ) 


9 x 10 8 


Low ionizing flux at the left grid edge 


<t> j ( rm — * c; — 1 "1 


4.5 x 10 9 


^Icdium ionizing flux at the left grid edge 


(f>i ■ f cm — ^ « — -M 
^hi b ) 


9 x 10 9 


ITicn irMi 171 n tr fliiY at tno 1 fift crnH pHtrp 

lllj^ll IMlllZjlllli 11UX dl 111L, Hjll til 1U ljUtlt_. 


L(pc 3 ) 


4.87 3 


Grid size 


^cclls 


128 3 


Number of grid cells 


^pix 


401 2 


Synthetic image pixels 


W) 


2.5 


Angular width per pixel 


N 1 


10 8 


Photons packets used in synthetic image 


D (pc) 


1000 


Observer distance 


dust / gas 


1 x 10~ 2 


Dust to gas ratio 


logi (He/H) 


-1 


Helium abundance 


log 10 (C/H) 


-3.66 


Carbon abundance 


logio(N/H) 


-4.40 


Nitrogen abundance 


logi (O/H) 


-3.48 


Oxygen abundance 


logic, (Ne/H) 


-4.30 


Neon abundance 


log 10 (S/H) 


-5.05 


Sulphur abundance 



ter the starting density distribution for this model, accu- 
mulating only a small amount of material before achieving 
pressure balance and establishing a weak photo-evaporative 
flow. As a result the enhanced cooling due to forbidden line 
radiation and the more complex thermal balance calcula- 
tion have increased the size of the neutral gas region, i.e. 
the Stromgren radius is smaller. The resulting imaged ob- 
ject has the curvature of a class A BRC, though it is not 
actually bright rimmed, beneath which hints of the old neu- 
tral cloud structure can be seen. 



5.3 Neutral gas properties 

Spectral energy distributions which include dust and free- 
free emissivities were generated for each converged grid, 
these are shown in Figure |10| Short of the near UV regime 
(from approximately 380 nm bluewards) scattered stellar 
photons dominate the signal and the flux decreases in ac- 
cordance with the star's distance from the grid, i.e. the 
strongest signal is from the high flux model and the weak- 
est is from the low flux model. Over most of the rest of 
the spectrum, where direct thermal radiation dominates the 
SED, the medium flux model consistently has the strongest 
signal, followed by the low then high flux models. At wave- 
lengths greater than about 370 (jm the low flux spectrum 
has the strongest signal due to the greater spatial extent 
of the cloud and therefore more widespread dust emission. 
The 10 ;Um silicate feature is clearly visible in emission in all 
SEDs. 

The dust temperatures and cloud masses are derived 
using the greybody fitting method described in section [3~T| 
We use a value of 214 gem -2 for C v and fit the spectrum 
between 450 and 850 /im following Thomps on et al.j ( |2004[ ). 
The resulting temperatures and cloud masses for each sys- 
tem are given in Table [5] The temperatures are all close to 
the dominant (that for at least 95% by mass) neutral cloud 
temperature of 10 K. The derived temperature is increased 



slightly by 1 — 2K by warmer dust at the interior edge of 
the bright rim. 

The calculated neutral mass in the low and medium 
flux cases are overestimated by a factor of approximately 2 
and 3.6 respectively. This is due to the standard practice 
of using a fixed value of C„ (e.g. Thompson et al. 2004 
Morgan et al.|2008 1 which is probably too large in this case 
The correct mass would have been inferred using values of 
110 and 60 gem -2 for C v in the low and medium flux mass 
calculations respectively. The inferred mass of the high flux 
cloud is close to the known value on the computational grid. 
This is because lower densities were attained in this model 
making the adopted value of C v more appropriate. 

Cloud masses calculated u sing the commonly ado pted 
value of 50 gem" 2 for C„ (e.g. Thompson et al. 20041 are 



also given in Table [5] These give a reasonable estimate of 
the medium flux mass but underestimate the low flux mass 
by over a factor of two. The large error in the high flux mass 
when using the lower value of C v is less significant due to 
the differences discussed in section [5~T1 

The sensitivity of this calculation to the dust temper- 
ature at low values is also illustrated in Table [5j where the 
calculated mass assuming a cloud temperature of 10 K (the 



dominant temperature in the neutral cloud, see section 5.1 1 



is also given for both values of C v used here. In the low and 
high flux cases the 1 K variation in temperature modifies the 
calculated masses by a factor in the range 8 — 16 %. The 2 K 
variation in the medium flux case leads to a change in the 
calculated mass by about 35 %. 

This mass calculation is already treated with caution 



(e.g.|Kerton et al.|200"T||Thompson et al.|2004||Morgan et al 



2008 1 these results suggest that using a fixed value of C„ for 



a range of clouds of different class will induce errors of a 
factor up to around 3.6 in some of the cloud masses and 
that the calculation is sensitive to temperature at values 
around 10 K with a typical difference in the mass of about 
15 ^oK" 1 . 
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Table 5. Neutral cloud properties from SED fitting. Included are two choices of C„ and the masses calculated using the dominant 
temperature of 10 K. 



Model 


Dust 


Cloud 


mass, (Mq) 


Cloud mass, (M©) 


Cloud mass, (Mq) 


Cloud mass 




temperature (K) 


C v = 


214gcm -2 


C„ = 50 gem -2 


C„ = 214/50gcm" 2 , T = 10K 


from grid (Mq) 


Low flux 


11 




144 


34 


159/37 


74 


Medium flux 


12 




75 


18 


115/27 


21 


High flux 


11 




76 


18 


90/21 


77 



Table 6. HPBW sizes and pixel noise levels for the images based 
on 30 s VLA exposures. 



Configuration 


HPBW 


rms pixel noise (mjy pixel x ) 


B 


3.9" 


0.27 


C 


12.5" 


2.7 X 10~ 2 


D 


44" 


2.2 X 10" 3 



5.4 Radio analysis 

20 cm radio images are used to determine the photoioniz- 
ing flux, IBL electron density and mass loss rate following 
the discussion in section [372] Synthetic images at 20 cm are 
generated, smoothed and subjected to noise in the man- 
ner described in sections |2.3| and |2.4| to be representative 
of the resolution and noise level of 30 second exposures of 
the VLA type B, C and D configurations. 30 second expo- 
sures using the D configuration are typical of those used in 
the NRAO VLA Sky Survey (NVSS) ( [Condon et aL]p98l ). 
We assume a bandwidth of 43 MHz and the use of 26 an- 
tennae. The rms noise associated with these observations 
is 0.524 mjybeam -1 . Table [6] summarises the HPBW value 
and rms pixel noise for each configuration. The regions used 
to derive the integrated IBL flux are the circular regions 
labelled 'B' on Ha images in Figure [12"| These all have an 
angular diameter of 1'. The cloud radii are those of the cir- 
cular regions marked 'A' on Figure [L2| and are 1.5, 0.6 and 
1.3 pc for the low, medium and high flux models respectively. 

The raw radio images, as well as simulated B, C and 
D configuration, 30 second VLA images are given in Fig- 
ure [TT] The colour scale used in each frame is set to match 
that of the raw image. As the beam size increases, the BRC 
structure is smoothed out and the brightness of the object is 
also modified. The calculated ionizing fluxes, electron den- 
sities and mass loss rates based on each radio image are 
summarised in Table [7] For each model, using a larger beam 
reduces the measured flux resulting in underestimates of the 
cloud properties by up to 25% relative to the unsmoothed 
image. This is because the compact, bright object flux is 
partially lost to the surroundings, which are uniformly much 
dimmer than the BRC. 

The ionizing fluxes at the left hand edge of the compu- 
tational grid are known to be, from low to high flux, 9 x 10 s , 

4.5 x 10 9 and 9 x 10 9 cm s . The tabulated values of ion- 
izing flux at the BRC, which are calculated at distance of 
around 2pc from the left hand edge of the grid, are there- 
fore of realistic magnitude. There is a discrepancy in that 
the high flux model ionizing fluxes are only slightly larger 
than the medium flux, however this is due to the change in 
ionization structure mentioned in section [FIT] that gives rise 



to an absence of a clear IBL and hence a lower measured 
flux. 

The electron density in the medium flux case has been 
correctly inferred as being significantly higher than that in 
the low flux case. The actual number densities in the IBL 
are in the range 60 - 190, 200 - 1500 and 50 - 100 cm" 3 
for the low, medium and high flux models respectively. In 
the medium and low flux cases the calculated values lie be- 
neath this range, suggesting that the effect of the IBL on 
compression might be underestimated. This underestimate 
arises due to cooler gas in the regions over which the flux 
is being integrated at the interior of the bright rim. These 
relatively cool regions arise where there is localised shielding 
from the stellar radiation field and radiative heating occurs 
primarily from the diffuse field. This is the cause of the rib- 
bing effect seen in the medium flux images, where the bright 
contours are those cells directly exposed to the stellar radi- 
ation field. Increasing the region over which the flux is inte- 
grated amplifies this problem. For example, in the extreme 
case of integrating the flux over the entire cloud (region A in 
Figure 12 l the low flux electron density is reduced to around 



15 cm . Conversely, reducing the size over which the flux 
is integrated too much will make the derived values more 
susceptible to noise. 

The inferred mass loss rates are higher for the larger 
clouds, which have a larger surface over which material can 
be lost. Providing an actual mass loss rate from the com- 
putational grids for comparison is non-trivial, particularly 
in the low flux case where the boundary between the BRC 
cloud and the gas in the wings of the model is poorly de- 
fined. We therefore calculate the difference between the mass 
at 195 and 200 kyr in a number of volumes encapsulating the 
BRCs on the computational grid to provide a range of mass 
loss estimates. These estimates are also shown in Table [7j 
the high flux model is neglected due to the changes to the 
ionization structure discussed in section 15.11 Both the low 
and medium flux mass loss estimates from the grid span just 
over one order magnitude, encompassing the mass loss rates 
from the 20 cm emission analysis. 

Although in the correct (but broad) range, the mass loss 
rate as a quantity does not give a proper indication of the rel- 
ative strengths of the photo-evaporative flows. For example, 
the medium flux model ejections are clearly more energetic 
than the low flux ones, giving rise to stronger rocket-motion, 
more rapid compression of the cloud and leaving a stronger 
signature in the ambient HII region, but its mass loss rate 
is similar. A more useful parameter for studying relative 
photo-evaporative flow strengths is the mass flux 



M t = 



1 M 



(15) 



n r 2 

where R is the radius of the cloud and f2 is the solid angle on 



Triggered star formation diagnostics 1 1 




Figure 7. Composite colour images of the low, medium and high 
flux models from top to bottom using Ha (red), the 5007 A [O 
III] line (green), and the 3968 A [Ne III] line (blue). Each image 
side spans 4.87 pc. 




Figure 8. A two colour image of the low flux model for an ob- 
server inclined thirty degrees relative to the planar ionization 
front. This image is constructed using the electron density sensi- 
tive 3729 A (red) and 3726 A (green), [O II] lines. The image side 
spans 4.87 pc. 




Figure 9. A two colour image of the medium flux model for an 
observer inclined thirty degrees relative to the planar ionization 
front. This image is constructed using the electron density sensi- 
tive 3729 A (red) and 3726 A (green) [O II] lines. The image side 
spans 4.87 pc 



the cloud bounded by the IBL. For near edge-on inclinations 
Q can be approximated assuming cylindrical symmetry. Ex- 
ample mass fluxes are also included in Table [7] which (based 
on the opening angle of the bow) assume that the low flux 
IBL bounds 2n(l — cos(n/3)) of the cloud surface and the 
medium flux cloud IBL bounds 2n of the surface. These 
mass fluxes imply that the medium flux flow is significantly 
stronger than that of the low flux model, in agreement with 
the behaviour from the model grid. The difficulty in using 
the mass flux as an indicator of relative photo-evaporative 
flow strengths is determining what fraction of the cloud is 
covered by the IBL. A first order comparison of the relative 
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Figure 10. A a comparison of the low, medium and high flux 
model spectra. 



mass fluxes of clouds of similar type could be obtained by 
assuming that the fraction of the cloud bounded by the IBL 



is the same between types (Sugitani et al. 19911 and that 



only the radius of the cloud varies. 



5.5 Diagnostic line ratio analysis 



As discussed in section [3T4] we emulate slit spectroscopy by 
calculating images at specific wavelengths and derive the 
conditions in the slit region using diagnostic line ratios. The 
slits are shown in the upper left of Ha images of the models 



in Figure 12 (marked C). Each slit has a width of 2 pixels 
and a length of 100, 50 and 80 pixels for the low, medium 
and high flux models respectively. These are placed over the 
ionized cells closest to the tip of the BRC. 

The electron densities and the average of the tempera- 
ture diagnostics are given in Table [8] along with the values 
calculated using radio diagnostics and those from the model 
grids. We retain the level of accuracy determined in sec- 
tion [471] As with the radio diagnostic, the electron density 
is clearly underestimated in the medium and low flux cases 
(e.g. hydrogen number densities are of order 60 — 190 cm -3 
for the low flux model). These values could be influenced 
by the low density (of order ten hydrogen atoms per cu- 
bic centimetre) region just beyond the BRC that has been 
excavated by the photo-evaporative flow. 

The temperatures derived using the [O III], [N II] and 
[Ne III] ratios are also given in Table [8] The temperature 
calculations use the electron densities calculated from the 
[O II] lines, but depend only weakly on this value so the 
largest error in n e (An e = 1460 in the medium flux case) 
translates into an error of only 15 K. The separate temper- 
ature diagnostics of each model are within 110K of one an- 
other and all lie within the range seen in the temperature 
maps of Figure [6] The actual temperature in the IBL is typ- 
ically 8000-9000 K so the average temperatures in Table [8] 
are more accurate than assuming a value to 10 4 K. Given 
that the electron densities are underestimated by both ra- 
dio and diagnostic line ratio techniques and the inferred IBL 
temperatures are more accurate using diagnostic line ratios, 
this technique is a viable tool for the study of RDI. 




Figure 12. Ha images of each model, low to high flux from top 
to bottom. Overlaid are the regions which are designated as the 
cloud cross section (those labelled A), the regions over which the 
radio flux is integrated (those labelled B) and the spectroscopic 
slits (those labelled C). 



5.6 Identifying RDI 

If the BRCs are in pressure equilibrium with the IBL then 
the virial theorem gives 



P 



3c: M c 



3GA/ C 



(16) 



where M c , R c , c s and Pi are the cloud mass, cloud radius, 
cloud gas sound speed and pressure in the IBL respectively 
(e.g. |Hartmannp5o"9"] |Thompson et al. 12004] ). If the BRC is 
not in pressure equilibrium then the contracting or expand- 
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Figure 11. Radio images of the low, medium and high flux models from top to bottom. The left hand panels are the raw output from 
TORUS. The other panels, from left to right, are the raw image modified to be respresentative of 30 s exposures using the VLA type B, C 
and D configurations. The grey scale of each image is set to match that of the raw output image. Each image side spans 4.87 pc 



Table 7. BRC IBL properties from analysis of radio images at 20 cm. The mass loss rate from the model grid is also included. 



Model 


Configuration 


Ionizing flux 


Electron density 


Mass loss rate 


Model mass los 


s rate 


Mass flux 






(<!>, cm -2 s _1 ) 


(n c , cm -3 ) 


(M, Mokyr" 1 ) 


(M, Mq kyr 


- 1 ) 


(M f ,M© Myr-ipc" 2 ) 


Low flux 


raw 


5.0 X 10 8 


45 


0.18 


1.4 x 10~ 2 - 


0.2 


26 


Low flux 


B 


4.9 x 10 8 


15 


0.18 


1.4 x 10" 2 - 


0.2 


25 


Low flux 


C 


4.6 x 10 8 


43 


0.17 


1.4 x 10" 2 - 


0.2 


24 


Low flux 


D 


3.3 x 10 8 


36 


0.15 


1.4 x 10~ 2 - 


0.2 


21 


Medium flux 


raw 


1.5 x 10 9 


122 


7.8 x 10~ 2 


9.8 x 10"3 - 


0.1 


35 


Medium flux 


B 


1.4 x 10 9 


119 


7.7 x 10" 2 


9.8 x 10" 3 - 


0.1 


34 


Medium flux 


C 


1.3 x 10 9 


113 


7.3 x 10" 2 


9.8 x 10" 3 - 


0.1 


32 


Medium flux 


D 


8.2 x 10 8 


91 


5.9 x 10~ 2 


9.8 x 10" 3 - 


0.1 


26 


High flux 


raw 


1.7 x 10 9 


89 


0.27 








High flux 


B 


1.7 x 10 9 


89 


0.27 








High flux 


C 


1.6 x 10 9 


87 


0.26 








High flux 


D 


1.4 x 10 9 


80 


0.24 









ing nature of the cloud can be inferred from the relative 
magnitudes of the left and right hand sides of equation [l6j 
which will be referred to as the external and supporting 
pressures (Pi and P s ) respectively. 

For the radio method we use the electron densities 
calculated in section |5.4| and make the standard assump- 



tion that the ionized gas is at 10000 K. For the diagnos- 
tic ratio method we use the electron densities and the av- 
erage temperatures calculated in section |5.5| The neutral 
cloud properties are those determined in section |5.3| using 
C„=214gcm~ 2 . Under this virial equilibrium technique a 
cloud radius needs to be assumed, the circular regions on 
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Table 8. Diagnostic line ratio electron densities and average IBL temperatures. The average electron densities calculated using the radio 
diagnostic are also included as well as the known conditions from the model. 



Model 


Line ratio n e (cm 3 ) 


Radio n a 


(cm 3 ) 


Model n c 


(cm 3 ) 


Line ratio T c (K) 


Model T c (K) 


Low flux 


20 


36 - 


45 


60 - 


190 


8500 


8000 - 9000 


Medium flux 


40 


91 - 


122 


200 - 


1500 


8500 


8000 - 9000 


High flux 


60 


80 - 


91 


50 - 


100 


8380 


8000 - 9000 



Figure [12] represent the cross section of the spherical clouds 
used to investigate equilibrium. 

The low, medium and high flux clouds have radii of 1.5, 
0.6 and 1.3 pc respectively. Given that most of the neutral 
mass on the grid is concentrated in the BRC, we assume 
in each stability analysis that all of the SED fitted mass 
is contained within the highlighted regions. The support- 
ing cloud pressures for each model, from low to high flux 
are 4.6 x 1CT 13 , 4.1 x 1CT 12 and 3.7 x 10 _13 dyncm~ 2 re- 
spectively. Typical pressures in the neutral cloud on the ra- 



diation hydrodynamic simulation grid are 10 



10" 



10" 12 - 10~ 10 and 10~ 13 - 10~ 12 dyncm- 2 in the low, 
medium and high flux models respectively so the derived 
pressure values are towards the centre of this range in the 
low and high flux models and towards the lower end of the 
range found in the medium flux model. 

The external pressures and their ratio relative to the 
support pressure are given in Table [9] The pressures in the 
IBL from the model grids range from around 0.5 — 3 x 10 -10 , 
l-5xl0~ 10 and 1-2. 5x 10" 10 dyncm" 2 for the low, medium 
and high flux models respectively. The IBL pressures derived 
using diagnostic line ratios lie within this range, however be- 
cause the assumed temperature of 10 4 K is larger than the 
actual temperature in the IBL the radio diagnostics overes- 
timate the pressure in the IBL in the medium and high flux 
cases. The low flux electron density was sufficiently underes- 
timated for the pressure to fall into the range on the model 
grid. 

All diagnostics imply that the cloud should be undergo- 
ing thermal compression, which is in qualitative agreement 
with the known model behaviour in the medium and low 
flux cases. Furthermore, the difference in pressure is weaker 
in the medium flux case than the low flux case. This reflects 
the relative behaviours of the clouds, as the rate of com- 
pression of the medium flux cloud has started to plateau 
by this point (200 kyr) in the radiation hydrodynamics cal- 
culation. The biggest pressure differences are found in the 
high flux model. Although this is not in agreement with the 
pressure supported system at the time of imaging in the ra- 
diation hydrodynamics calculation, the ionization front has 
been relocated to a point at which D-type driving of the 
front into the cloud is about to occur as discussed in section 
|5.1| Given these modifications to the starting conditions, it 
is likely that the resulting BRC following a new calculation 
with metals would more closely resemble that of the medium 
flux model. The ratio of external to supporting pressure lies 
within the range of those found on the computational grid 
in all cases, but are typically towards the lower end of the 
range. This is due to a combination of the overestimated 



neutral cloud masses and the underestimated IBL electron 
densities. 



6 SUMMARY AND CONCLUSIONS 

We have performed a range of synthetic observations and 
diagnostics of radiation hydrodynamic RDI model results 
to investigate the accuracy of diagnostic techniques and to 
identify signatures of RDI. We have produced SEDs, images 
that are representative of 30 second VLA radio observations 
using the type B, C and D configurations (the latter of which 
is used for the NRAO VLA Sky Survey) and forbidden line 
images of each model. These synthetic observations have 
been used to replicate a number of diagnostics to calculate 
the conditions in the neutral BRC and its IBL. Using these 
conditions we have performed a virial stability analysis of 
each system to determine whether the BRCs are being com- 
pressed by the IBL. 

We draw the following conclusions from our imaging 
and diagnostics: 

1. The synthetic images generated show objects which 
are morphologically similar to BRCs observed in star 
forming regions. 

2. The neutral cloud dust temperatures derived using 
greybody fitting of the SED are similar to those calculated 
in observational studies (e.g. |Thompson et al . 2004 ; Morgan 
et al. 2008 1 and are slightly higher than, but within 2 K 



of, the known dominant cloud temperature (10 K for at 
least 95% of the cloud by mass). This slight temperature 
overestimate is due to warmer dust in the cloud. The 



mass inferred using equation 10a from Hildebrand ( 1983 1 
(equation [6] in this paper) at 10 K is found to vary by 
8 — 16 % with a change of 1 K and by up to 35 % with 
a change of 2 K. Using a fixed value of C„ for clouds of 
different class is found to induce an error in the calculated 
mass of up to a factor 3.6. 

3. The temperatures derived using diagnostic line 
ratios at around 8500 K are more accurate than assuming a 
canonical temperature of 10 4 K as in the radio diagnostics. 
The electron number densities established using both 
techniques are underestimates of those in the IBL on the 
computational grid. Diagnostic line ratios underestimate 
the electron density due to contamination from the low 
density region excavated by the photo-evaporative flow. 
The radio electron density underestimate is due to lower 
emission from the interior regions of the BRC where 
some shielding from the stellar radiation field occurs and 
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Table 9. IBL pressures and their ratio to the supporting cloud pressure in virial stability analysis. The ratio of IBL to neutral cloud 
supporting pressures from the model are also included. 



Model 


Diagnostic method 


External pressure (Pi, xlO 10 dyncm 2 ) 


Pi/P s 


Model Pi/P s 


Low flux 


Radio, raw image 


2.1 


460 


50-3000 


Low flux 


Radio, VLA configuration B 


2.1 


460 


50 - 3000 


Low flux 


Radio, VLA configuration C 


2.0 


430 


50 - 3000 


Low flux 


Radio, VLA configuration D 


1.7 


370 


50 - 3000 


Low flux 


Slit spectroscopy 


0.78 


170 


50 - 3000 


Medium flux 


Radio, raw image 


5.6 


140 


1 - 500 


Medium flux 


Radio, VLA configuration B 


5.5 


130 


1 - 500 


Medium flux 


Radio, VLA configuration C 


5.2 


130 


1 - 500 


Medium flux 


Radio, VLA configuration D 


4.2 


100 


1 - 500 


Medium flux 


Slit spectroscopy 


1.6 


40 


1 - 500 


High flux 


Radio, raw image 


1.1 


1110 


100 - 2500 


High flux 


Radio, VLA configuration B 


4.1 


1110 


100 - 2500 


High flux 


Radio, VLA configuration C 


4.0 


1080 


100 - 2500 


High flux 


Radio, VLA configuration D 


3.7 


1000 


100 - 2500 


High flux 


Slit spectroscopy 


2.3 


620 


100 - 2500 



radiative heating is due to the, relatively weak, diffuse field. 

4. The supporting cloud pressures are found to lie 
within the range on the computational grids for all models. 
The IBL pressures are found to be within the correct range 
for the forbidden line ratio diagnostic. However, since the 
assumed temperature in the IBL is larger than the actual 
temperature the radio diagnostics overestimate the IBL 
pressure in the medium and high flux cases. The IBL 
pressure is greater than that of the supporting cloud for 
each model, implying that the cloud is being compressed. 
This behaviour is qualitatively correct for the low and 
medium flux models, however is not the observed behaviour 
of the pressure balanced high flux model. The ratio of 
external to supporting pressure is found to be towards the 
lower end of the range of values found on the model grid in 
each case. 

5. The ionization state of the low and medium flux 
grids following the initial photoionization calculation 
remained reasonably consistent with that at the end of 
the radiation hydrodynamic calculations. However, for 
the high flux model in which the density structure in 
the system was not significantly altered in the radiation 
hydrodynamics calculation, the ionization front has moved 
closer to the star. This is due to the enhanced cooling by 
forbidden line emission and more comprehensive thermal 
balance calculation, essentially resetting the model to a 
point in which D-type expansion of the ionization front 
is yet to occur. With a greater distance over which to 
accumulate material, the modified high flux system may 
well have achieved a stronger photo-evaporative flow and 
higher levels of compression of the BRC than in the original 
model. This result indicates the importance of including 
atomic chemistry in future calculations. 

6. The effect of moving to larger beam sizes in radio 
observations is to contaminate the integrated flux from 
the BRC with that from the surrounding neutral cloud 
and HII region. Since both the neutral cloud and HII 
region are dimmer than the IBL this reduces the integrated 



flux, resulting in an underestimate of the incident ionizing 
flux, electron density and mass loss rate relative to that 
calculated using an unsmoothed comparison image 

7. Despite the presence of strong photo-evaporative 
flows which establish prominent, relatively cool, low density 
regions on the model grid, they are difficult to detect 
in the images calculated here. No striations are detected 
in Ha about the BRC, rather the visual indicator of 
photo-evaporative flow in these systems is darker regions in 
the vicinity of the BRC, where material is excavated by the 
flow. It is likely that striations are not observed because 
the resolution in these calculations was not sufficiently high 
to resolve them. 

8. The estimate for mass loss rate of clouds from 



Lcfloch & Lazareff (1994) is found to be in good agreement 
with that of the clouds on the computational grid. By using 
the mass flux, a measure of the relative strengths of the 
photo-evaporative flows can be calculated, correctly imply- 
ing that the medium flux model exhibits much stronger 
flow than the low flux model, despite having a similar mass 
loss rate. 



We have shown that existing diagnostics do give an in- 
sight into the conditions of BRCs and can broadly be used 
to infer whether or not RDI is occuring using virial stability 
analysis. However, by comparing the inferred conditions to 
those on the model grids we find that these diagnostics are 
each subject to significant sources of error. The cumulative 
effect of these errors is IBL-to-cloud pressure ratios towards 
the lower end of the known range, suggesting that the effects 
of shock compression might be underestimated. 

We next intend to perform synthetic diagnostics using 
molecular lines, which are a useful tool for investigating the 



kinematic behaviour of astrophysical gases (e.g. De Vries 



et al.| 


2002| 


Urquhart et al.|2006||Morgan et al.|2009||Rundle 


et al. 


2010 


. Applying these techniques we can draw com- 



parison conculsions to those in this paper from a more so- 
phisticated analysis of the gas, investigate the kinematic be- 
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haviour of BRCs and provide a test of widely used molecular 
line diagnostics. 
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